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Abstract : The study of convective heat and mass transfer flow 
over a vertical plate with nth order chemical reaction in a 
porous medium is analyzed both analytically and numerically. 
The resulting governing boundary layer equations are highly 
non-linear and coupled form of partial differential equations 
and have been solved by using implicit finite difference method 
of Crank-Nicolson. To check the accuracy of the numerical 
solution, steady-state solutions for velocity, temperature and 
concentration profiles are obtained by using perturbation 
method. The effects of various parameters such as thermal 
Grashof number, solutal Grashof number, chemical reaction, 
order of chemical reaction, radiation, thermal conductivity, 
Prandtl number and Schmidt number on the velocity, 
temperature and concentration profiles as well as skin friction, 
Nusselt number and Sherwood number are presented 
graphically. The numerical results showed that both 
temperature and concentration increases with the increase in 
chemical reaction while velocity decreases. It was also 
observed that the numerical and analytical solutions are found 
to be in an excellent agreement. 

Key words: heat and mass transfer, free-convection, porous 
medium and chemical reaction. 

1. Introduction 

The study of convective flow of heat and mass transfer with the 
influence of chemical reaction is given primary importance in 
science and engineering fields. This phenomenon plays an 
important role in chemical industry, hydro metallurgical 
industries, petroleum industry, cooling of nuclear reactors and 
packed -bed catalyst reactors Gregory and Ilhan [8]. 

Chemical reaction usually accompanies a large amount of 
endothermic and exothermic reactions. These characteristics 
could easily be seen in a lot of industrial processes. For example, 
in the power industry among the methods of generating electric 
power in which electrical energy is extracted directly from a 
moving conducting fluid. Recently, it has been realized that it is 



not always permissible to neglect the convection effects in 
porous constructed chemical reactors. The reaction produced in a 
porous medium was extraordinarily in common, such as in fuel 
cell modules and the polluted underground water because of the 
toxic substance Nield and Bejan [19]. 

In a recent literature review, porous media are widely used in 
high temperature heat exchangers, turbine blades, jet nozzles, 
solar collectors, drying processes, building constructions etc. In 
practice, cooling of porous structure is achieved by forcing the 
liquid or gas through capillaries of solid. Actually, they are used 
to insulate a heated body to maintain it is temperature. The effect 
of variable permeability on combined free and forced convection 
in porous media was analyzed by Chandrasekhar and 
Namboodiri [5]. Combined heat and mass transfer problems with 
chemical reaction are also important for a variety of engineering 
applications. These applications include solar collectors, grain 
storage and drying processes, heat exchangers, geothermal and 
oil recovery, building construction, nuclear waste material and 
many others Minto et al. [17]. In view of it is applications 
several authors have studied nth order chemical reaction and 
porous medium both in mechanical and physiological situations. 
Ahmad et al. [1] presented effect of first order chemical reaction 
and radiation on an unsteady MHD flow of an incompressible 
viscous electrically conducting fluid with variable temperature 
and mass transfer. Krishnendu [13] analyzed a steady boundary 
layer slip flow and mass transfer with nth order chemical 
reaction past a porous plate embedded in a porous medium. Also 
work on steady mixed convection flow in an incompressible 
viscous electrically conducting fluid with nth order chemical 
reaction has been reported by Gurminder et al. [9]. Additionally, 
a numerical investigation on the effect of chemical reaction on 
unsteady natural convection boundary layer flow over a semi- 
infinite vertical cylinder is carried out by Kawala and Odda [11]. 

A finite difference solution of the homogeneous first order 
chemical reaction on unsteady flow past an impulsively started 
semi-infinite vertical plate with variable temperature and mass 
transfer in the presence of thermal radiation have been examined 
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by Loganathan [14]. Mansour et al. [16] described the influence 
of chemical reaction and viscous dissipation on MHD natural 
convection flow. Recently, Dilal and Sewil [6] have studied a 
numerical model of the analysis of combined effects of mixed 
convection MHD heat and mass transfer in a visco-elastic fluid 
in a porous medium with chemical reaction, non uniform heat 
source/sink and visco-ohmic dissipation. Most recently, Saleh et 
al. [23] analyzed convective heat and mass transfer 
characteristics of an incompressible MHD visco-elastic fluid 
flow immersed in a porous medium over a stretching sheet with 
chemical reaction and thermal stratification. 

In another article, Al-Rashdan et al. [2] presented an analytical 
investigation to the problem of fully developed natural 
convective heat and mass transfer through porous medium in a 
vertical channel in the presence of a first order chemical 
reaction. El-Sayed [7] has presented the theoretical study of the 
effect of mass diffusion of chemical species with first order 
reaction on peristaltic motion of an incompressible Jeffery fluid. 
Additionally, the effect of chemical reaction on the forced and 
free convection boundary layer that flows in a semi-infinite 
expanse of an electrically conducting viscous incompressible 
fluid past a vertical porous plate is carried out by Rajeshwari et 
al. [22]. 

Analysis of transport processes and their interaction with 
chemical reaction has the greatest contributions to many 
chemical sciences. The chemical reaction effect on heat and 
mass transfer flow along a semi-infinite horizontal plate has 
been investigated by Anjalidevi and Kandasamy [3]. 
Additionally, Patil and Ali [21] analyzed heat and mass transfer 
from mixed convection flow of polar fluid along a plate in 
porous media with chemical reaction. A detailed numerical 
study has been carried out for unsteady hydromagnetic natural 
convection heat and mass transfer with chemical reaction over a 
vertical plate in rotating system with periodic suction by Parida 
et al. [20]. Mahdy [15] have considered effect of chemical 
reaction and heat generation or absorption on double diffusive 
convection from vertical truncated cone in a porous media with 
variable viscosity. Recently, a numerical investigation on the 
effects of chemical reaction, radiation and magnetic field on the 
unsteady free convection flow, heat and mass transfer 
characteristics in a viscous, incompressible and electrically 
conducting fluid past an exponentially accelerated vertical plate 
is carried out by Kishore et al. [12]. Most recently, Karunakar et 
al. [10] investigated effects of heat and mass transfer on MHD 
mixed convection flow of a vertical surface with radiation, heat 
source/absorption and chemical reaction. They obtained the 
results using perturbation method. 
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This work is an extension of Muthucumaraswamy and 
Manivannan [18] which studied unsteady free convective flow 
of a viscous incompressible flow past an infinite isothermal 
vertical oscillating plate with variable mass diffusion taking into 
account the homogeneous chemical reaction of first order. The 
mathematical formulations have taken into cognizance magnetic 
field, permeability in the momentum equation, variable thermal 
conductivity and radiation in the energy equation which are 
absent in Muthucumaraswamy and Manivannan [18]. Also, the 
boundary conditions for the flow are different. 

Hence, in this present study, special attention is given to 
convective heat and mass transfer flow over a vertical plate with 
nth order chemical reaction in a porous medium. 

2. Mathematical formulation 

Consider a free convection unsteady flow of an incompressible, 
viscous, electrically conducting and chemical reacting fluid with 
radiative heat and thermal conductivity past a vertical plate 

through a porous medium. The x - axis is taken to be along the 

plate and the y - axis is taken normal to the plate. A transverse 

magnetic field of strength B is applied normal to the y-axis. It 

is assumed that the induced magnetic field, the external electric 
field and the electric field are negligible. This assumption is 
justified by the fact that the magnetic Reynolds number is very 
small. The fluid properties are assumed to be constant. It is also 
assumed that there exists a homogeneous chemical reaction 
between the fluid and species concentration. The concentration 
of diffusing species is very small in comparison to other 

chemical species, the concentration of species from the wall, C 

is infinitely small and hence the Soret and Dufour effects are 
neglected Sarada and Shankar [24]. 

Under these assumptions and along with Boussinesq's 
approximation, the governing equation describing the 
conservation of mass, momentum, energy and concentration are 
described by the following system of equations respectively: 

8v 
oy 

du du ' d 2 u ( crB 2 v | 

V o = V ~ - H 7 U 

8t' 8y' dy 2 y p k" J (2) 
+gj3(T'-T' ) + gj3 t (C-C ) 
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dC dC d C 

— -v B — = D — -R {C -CJ 

dt dy dy 



(3) 



dT' dT' _ k d 
dt ° dy pC p dy 
1 da 



\ + a(T\ r T\) 



dT 

dy' 



(4) 



pC„ dy 



The corresponding initial and boundary conditions are 
prescribed as follows: 



t < 0, u =0,T -> T , C -> C for all y 



t>0, u =0,T 



C = C at y = 



(5) 



0, T = T U ,C =C U at y 



From continuity equation, it is clear that the suction velocity is 
either a constant or a function of time. Hence, on integrating 
equation (1), the suction velocity normal to the plate is assumed 
in the form, 



where v is a scale of suction velocity which is non-zero positive 

constant. The negative sign indicates that the suction is towards 

the plate and v > corresponds to steady suction velocity 

normal at the surface. The third and fourth terms on the right 
hand side of equation (2) denote the thermal and concentration 
buoyancy effects respectively, The last term of equation (3) is 
the chemical reaction term and n is the order of chemical 
reaction which relates the rate of chemical reaction with the 
concentrations of the reacting substances and the higher the 
order of the reaction the smaller is the region in the parameters 
for which multiplicity can occur, while the last term of equation 
(4) represents the radiative heat flux term. 

The Radiative heat flux q under Rosseland approximation by 
Brewster [4] has the form: 



4<y df 

3k dy 



(6) 



where <t is Stefan- Boltzmann constant and k is the mean 
absorption coefficient. 

We assume that the temperature differences within the flow 
are so small that T 4 can be expressed as a linear function of the 
temperature. This is accomplished by expanding T 4 in a Taylor 
series about T Q and neglecting the higher order terms. Thus 

T 4 = 4T Y - 3T 4 (7) 
Using equations (6) and (7) in equation (4) we obtain: 



dt' 



dT' 
~dy 



a 



pC„ dy 



dT 

\ + a(T\-T' ) — 

dy 



16(77/ d'T 
3kpC dy' 2 



(8) 



where u and v are the Darcian velocity components in the x- 
and y- directions respectively, t is the time, v is the kinematic 
viscosity, g is the acceleration due to gravity, /? is the 
coefficient of volume expansion, p the density of the fluid, a is 

the scalar electrical conductivity, p is the volumetric 
coefficient of expansion with concentration, C is the specific 

heat capacity at constant pressure, k is the permeability of the 
porous medium, k is the dimensionless thermal conductivity of 
the ambient fluid, a is a constant depending on the nature of the 
fluid, R* is the dimensionless chemical reaction, n is the order 
of chemical reaction, D is the coefficient of molecular 
diffusivity, v is the constant suction parameter, q r is the 

radiative heat flux in the y- direction, B is the magnetic 

induction of constant strength. T and T are the temperature of 
the fluid inside the thermal boundary layer and the fluid 
temperature in the free stream respectively, while C and C 
are the corresponding concentrations. 

To obtain the solutions of equations (2), (3) and (8) subject to 
the conditions (5) in non-dimensional forms, we introduce the 
following non-dimensional quantities: 
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c c 

C = ° ,Pr = 



c c 



-^,5c=— . 



V" 



2 



-,Gr = 



(9) 



g ( C " C °) t \ 

Gc = ^ '-, = (T v T ), 

-_o ^ = ^L j? = 16a o C 



Applying these non-dimensionless quantities (9), the set of 
equations (2), (3), (8), and (5) reduces to the following: 



du du d u 

dt dy dy 2 



f 



1 

M + - 



k. 

+ Gr0 + GcC 

dc dc i d 2 c 

Y — = - KC 

dt dy Sc dy" 



(10) 



(11) 



d_6_ 
dt 



d_e__i_ 

dy Pr 
1 



Pr ay Pr 

With the following initial and boundary conditions 

t < u = 0, 6 = 0, C = 0, for all y, 
t > u = 0, G = 1, C = 1, at y = 0, 
M = 0, = 0, C = a? y = l. 



(12) 



(13) 



where Pr is the Prandtl number, 5c is the Schmidt number, M 
is the Magnetic field parameter, Gr is the thermal Grashof 
number, Gc is the solutal or mass Grashof number, X is the 
variable thermal conductivity, y is the variable suction 

parameter, R is the radiation parameter, K is the chemical 

reaction parameter, n is the order of chemical reaction, k is the 
permeability parameter, t is the dimensionless time while u 



and v are dimensionless velocity components in i - and y - 
directions respectively. 

3. Analytical solutions 

The governing equations involved in this problem are highly 
coupled and non linear. In order to confirm the validity of this 
model, it is therefore of interest to reduce the governing 
equations of this problem to a form that can be solved 
analytically. A special case of this problem that exhibits 
approximate solution is the problem of steady convective heat 
and mass transfer flow over a vertical plate with chemical 
reaction and porous media. The steady state equations and 
boundary conditions for the problem can be written as follows: 



d 2 u 



du 



r 



r- 



dy dy 



M 



V 



d 2 C 



kj 

Grd + GcC 

dC 



(14) 



2 - + ySc KScC 

dy dy 



d 2 d6 

— + /Pr R0 = O 

dy' dy 



The boundary conditions are: 



(15) 



(16) 



u = 0, = 1, C = 1, at y = 0, 
u = 0, = 0, C = at y = l. 



(17) 



To find the approximate solution to equations (14)-(16) subject 
to equation (17), we use perturbation method which is a method 
that is used to approximate the solution to a differential equation 
analytically. Therefore the physical variables u, C , and 6 can 

be expanded in the power of(7?«l). This can be possible 

physically as R for the flow is always less than unity. Hence we 
can assume solution of the form 

u = u (y) + Ru l (y) + 0(R 2 ) 

C = C (y) + RC l (y) + 0(R 2 ) (18) 

d = e (y)+Re l (y)+o(R 2 ) 

Using equation (18) in equations (14)-(17) and equating the 
coefficient of like powers of R, we have 
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v 



M +-\u =-Gr0 n -GcC n 



C +yScC -K r ScC =0 

o "+rPr0 o =0 



f 



u t + yu t 



1 



M + 

kj 



u. = -GrO, - GcC, 



C," + yScC l - K r ScC l = 

The corresponding boundary conditions are 

u a = 0, = 1, C = 1 

Ul = 0, (9j = 0, C, = a/ j = 

Mq =0, o =O, C =0 

M] = 0, (9j = 0, C, = at y = 1 



(19) 

(20) 
(21) 



(22) 
(23) 
(24) 



(25) 
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Solving equations (19)-(24) with the help of equation (25), we 
get 



u =L c /' y +L m e 



L l4 e 



L n +L l2 e 



-y Pry 



L l3 e 



(26) 



9 e +L 10 e- +L n +L 12 e ■ + L l3 e 



L l4 e' v +R 



J i- _i_ T 



-j- Pry 



C = L 7 e"' y + L^e" 2 " 



= L, +L,e rny +R 



L,+L,e~ rVty + 



(32) 



(33) 



(34) 



4. Numerical solutions 



The system of transformed coupled non-linear partial differential 
equations (10), (11) and (12) together with the initial and 
boundary conditions (13) have been solved numerically using a 
finite difference technique of implicit type namely Crank- 
Nicolson implicit finite difference method which is always 
convergent and stable. The equivalent finite difference 
approximations corresponding to equations (10)-(12) are given 
as follows: 



' u — u. . ^ 

I.J+l I, J 



At 



r 



2Ay 



2(Ay 2 ) 



u. — 2u . , + u . ., 

z+1,7+1 ',7 + 1 '—1,7+1 

, u . — 2u +u ,. 



M + - ]m w + Gr{e i . ) + Gc(C tJ ) (35) 



M l = V + +L 17 +L 18 e ^ PrV + 



(L 19 y + L 20 )e-^ 



a =o 



3, = A + V 



■/Pry 



(27) 

(28) 
(29) 
(30) 



9 l = L 3+ L 4 e-^ + (L 5 -L 6 )ye-^ (31) 
In view of the above equations the solutions are 



C . , -C 



2Ay 



(C + , -C , ) 



25c (Ay 2 ) 



^C,, . ,-2C. + 

1+1,7+1 i,7+l 
V 1-1,7+1 >+l,J >,7 1-1,7 / 



(36) 
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1,7+1 i, J 



At 
H 



r 



2Ay 



(ft, -9 , .) 



(37) 



'0. .. 20 , • 
2Px(Ay 2 ){9 +1J -29.+9 i _ lj 

—*^(e +l .-9..) 2 --(9.) 



The initial and boundary conditions take the following forms: 



u. . = 0, 0..=O, C =o 

«o.y=0, O , = 1, C .=l (38) 

« . = 0, (9, . = 0, C. .=o 

h,j h,j h,j 



where h corresponds to 1 . 

The index i refers to space y and j refers to time t. Ay and At are 

the mesh sizes along y - direction and time t - direction 
respectively. The Finite difference equations (35)-(37) at every 
internal nodal point on a particular n-level constitute a tri- 
diagonal system of equations which are solved by using the 
Thomas Algorithm. 

In each time step, the concentration and temperature profiles 
have been computed first from equations (36) and (37) and then 
the computed values are used to obtain the velocity profile at the 

end of time steps that is u. . +] computed from equation (34). This 

process is carried out until the steady state solution of the 
convergence criteria for stability of the scheme is assumed to 
have been reached. Computations are carried out for different 
values of physical parameters involved in the problem. 

The skin friction coefficient, the Nusselt number and Sherwood 
number at the plate are given by: 







^89} 








,Nu = - 




,Sh = 






y=0 




y=0 


K Q y J 



c. 



5. Results and discussion 



(39) 



In order to get a physical view of the present problem, numerical 
computations have been carried out for different values of 
thermal Grashof number ( Gr ), solutal Grashof number ( Gc ), 
Magnetic parameter ( M ), porous parameter ( k ), suction 



parameter ( y ), radiation parameter ( R ), chemical reaction 

parameter ( K r ), order of chemical reaction ( n ), variable 

thermal conductivity ( X ), Prandtl number ( Pr ) and Schmidt 
number ( Sc ). The purpose of this computation given here is to 
assess the effect of these physical parameters upon the nature of 
the flow and transport. Computations are obtained for fluids with 
Prandtl number (Pr = 0.71, 1.0 and 7.0) corresponding to air, salt 
water and water respectively. The diffusing chemical species of 
most common interest in air have Schmidt number and is taken 
for (Sc = 0.22 0.60) which corresponds to hydrogen and oxygen 
respectively. The default values of the physical parameters are as 
follows: 

Gr = 2.0, Gc = 2.0, M =1.0, K r = 2.0, Pr = 0.7 1, 
Sc = 0.22, R = l.0,k = 0.1,n = 0A,y = 0.2, A = 1.0. 

All graphs therefore correspond to these values unless 
specifically indicated on the appropriate graphs. 




Fig. 1. Steady- state solution of Velocity profile for 
numerical and analytical solutions. 
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Fig. 2. Steady-state solution of Temperature profile for 
numerical and analytical solutions. 
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1 Steady sIbIb nf Numerical solution 
— Steady stats nf Analytical solutSnri 




Fig. 3. Steady-state solution of Concentration profile for 
numerical and analytical solutions. 

Figures (l)-(3) present comparison of the steady-state solution 
for numerical and analytical solutions of velocity, temperature 
and concentration profiles respectively.These results are 
presented to illustrate the accuracy of the numerical and 
analytical solutions. It is observed that the time necessary to 

reach the steady-state solution is t m „ , where f v is the steady- 

J UldX max - 1 

state value of time (t). Hence, the numerical solution for the 
problem is in good agreement with analytical solution at time 



Figures (4) to (22) are all plotted at time t = 0.2 respectively. 




Fig. 4. Velocity profile for different values of 
thermal Grashof number ( Gr ). 




Fig. 5. Velocity profile for different values of solutal 
Grashof number ( Gc ). 

Figures (4) and (5) describe typical velocity profiles in the 
boundary layer for various values of thermal Grashof number 
and solutal Grashof number in cases of cooling of the surface 
respectively. It is observed that an increase in Gr and Gc 
results in an increasing momentum boundary layer, this is due to 
the fact that increasing the values of Gr and Gc has the 
tendency to increase the thermal and mass buoyancy effect. This 
gives a rise to an increase in the induced flow. 




Fig. 6. Velocity profile for different values of 
porosity parameter ( k ) 
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Fig. 7. Velocity profile for different values of 
Magnetic parameter ( M ). 

Figure (6) depicts the variation of velocity distribution across the 
boundary layer for various values of porous parameter k . The 
porosity medium has considerably effect on the velocity. An 
increase in porosity leads to increasing velocity profiles. The 
effect of magnetic parameter M on the velocity profile in the 
boundary layer is displayed in figure (7). Application of a 
transverse magnetic field produces a resistive-type force called 
the Lorentz force. This force has the tendency to slow down the 
motion of the fluid in the momentum boundary layer. 




Fig. 8. Velocity profile for different values of Prandtl 
number ( Pr ) 




Figures (8) and (9) describe the behavior of Prandtl number on 
velocity profile and temperature distribution. It is seen that 
momentum and thermal boundary layer both decrease with the 
increase in Prandtl number. This is because smaller values of Pr 
are equivalent to increasing the thermal conductivities and 
therefore heat is able to diffuse away from the heated plate more 
rapidly than for higher values of Pr . 
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Fig. 10. Velocity profile for different values of Radiation 
parameter ( R ) 




Fig. 9. Temperature profile for different values of Prandtl 
number ( Pr ) 



Fig. 11. Temperature profile for different values 
of Radiation parameter ( R ) 

The effect of radiation parameter R on the velocity and 
temperature profiles are plotted in figures (10) and (11). It is 
observed that as the value of radiation parameter increases the 
velocity and temperature profiles decreases with a decreasing in 
the momentum and thermal boundary layer thickness. 
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Fig. 12. Velocity profile for different values of 
Schmidt number ( Sc ) 
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Fig. 13. Concentration profile for different 
values of Schmidt number (Sc) 

In figures (12) and (13) the influence of Schmidt number Sc on 
velocity and concentration are displayed respectively. As the 
Schmidt number increases, the velocity and concentration 
decreases. This causes the concentration buoyancy effect to 
decrease yielding a reduction in the fluid velocity and 
concentration. Physically, increase of Sc means decrease in 
molecular diffusivity. 
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Fig. 14. Velocity profile for different values of 
Suction parameter ( y = q) 
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Fig. 15. Temperature profile for different values of 
Suction parameter ( y = q) 

Figures (14) and (15) represent respectively the velocity and 
temperature profiles for various values of suction parameter. It is 
noticed that the velocity and temperature profiles decreases 
monotonically with increasing suction parameter which 
indicating the usual fact that suction parameter stabilizes the 
boundary layer growth. 




Fig. 16. Velocity profile for different values of 
Thermal conductivity ( A, = A) 




Fig. 17. Temperature profile for different values of 
Thermal conductivity ( X = A) 
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Figures (16) and (17) illustrate the effect of thermal conductivity 
parameter on the velocity distribution and temperature profile 
respectively. Increase in the mixed convection parameter leads 
to increasing in the fluid flow and thermal boundary layer 
respectively. 
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Fig. 18. Velocity profile for different values of 
Chemical reaction parameter ( K r ) 
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Fig. 19. Temperature profile for different values of 
Chemical reaction parameter (K r ) 




Fig. 20. Concentration profile for different values of 
Chemical reaction parameter ( K ) 



Figures (18)-(20) depicts the influence of chemical reaction 
parameter on the velocity, temperature and concentration 
profiles in the boundary layer respectively. As shown the 

velocity and concentration are decreasing with increasing K , 

but the temperature profile increases as K increases, this is due 

to the fact that destructive chemical reduces the solutal boundary 
layer thickness. 




Fig. 21.Velocity profile for different values of order 
Chemical reaction ( n ) 




Fig. 22. Concentration profile for different values of 
order Chemical reaction ( n ) 

The effects of order chemical reaction ( n ) on the velocity and 
concentration profiles are displayed in figures (21) and (22) 
respectively. It is observed from these figures that an increase in 
the chemical reaction order increases the velocity and 
concentration profiles respectively. 
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Fig. 23. Skin friction for different values of t for 
Pr = 0.71 
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Fig. 24. Skin friction for different values of t for Pr =7.0 

In figures (23) and (24), skin friction is plotted against magnetic 
parameter M for different values of non-dimensional time ( t ) 
and Prandtl number ( Pr = 0.71,7.0 ). It is observed that as 
magnetic parameter and time increases the steady-state of the 
skin friction increases for both air (Pr=0.71) and water 
( Pr = 7.0 ) respectively. 




Fig. 26. Nusselt number for different values of t for 
Pr=7.0 

Figures (25) and (26) show representative profiles for Nusselt 
number (rate of heat transfer) for different values of non- 
dimensional time (?) and thermal conductivity ( A = A) for 
Prandtl number ( Pr = 0.71, 7.0 ). These figures reflects that an 
increase in time t and A increases the rate of heat transfer, this is 
physically true because as A increases the fluid temperature 
increases and consequently there is high rate of heat transfer on 
the boundary. 
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Fig. 25. Nusselt number for different values of t for 
Pr=0.71 



Fig. 27. Sherwood number for different values of t 
for Sc = 0.22 



-14r 

-16 



E 

P -18- 



1 -20- 



-22. 



-24- 



1.5 



2.5 
n 



— 1=10 
"—1=1.2 
■—1=1.4 
— SS 



3.5 



Fig. 28. Sherwood number for different values of t for 
5c = 0.60. 
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Figures (27) and (28) describe the behavior of Sherwood number 
against order of chemical reaction ( n ) with different values of 
time t for Schmidt number ( Sc — 0.22, 0.60). It is observed 
that as n and t increases the Sherwood number increases for 
hydrogen (Sc— 0.22) and decreases for water vapor 
(St =0.60). 
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Appendix 
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6. Conclusion 

The present numerical study has been carried out for convective 
heat and mass transfer flow over a vertical plate with nth order 
chemical reaction in a porous medium. An implicit finite 
difference method of Crank-Nicolson is employed to solve the 
equations governing the flow. The following conclusion has 
been drawn for the present numerical investigation. 

1. It is found out that the velocity profiles increases with 
increasing thermal Grashof number, solutal Grashof number, 
porous material, order of chemical reaction and thermal 
conductivity. It is also interesting to note that increase in 
magnetic parameter, chemical reaction, Prandtl number, 
radiation, Schmidt number and suction leads to decreasing the 
flow velocity. 

2. An increasing Prandtl number, radiation parameter, 
chemical reaction and suction decreases the temperature profiles 
of the fluid flow field at all points. However, the temperature 
increases with increasing thermal conductivity and chemical 
reaction. 

3. The concentration profile increases with increasing 
value in order of chemical reaction whereas it decreases with the 
increasing values of Schmidt number and chemical reaction. 

4. Skin friction increases with increasing magnetic field 
and time t. 

5. The rate of heat transfer in terms of Nusselt number 
increases with increasing thermal conductivity and time. 

6. An increasing order of chemical reaction increases 
Sherwood number for hydrogen (Sc = 0.22) and decreases for 
water vapor (Sc = 0.60). 

7. The numerical solution is in excellent agreement with 
the analytical solution at time . 
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Nomenclature 

C - Concentration 

C - Specific heat at constant pressure 

D - Mass diffusivity 

g - Acceleration due to gravity 
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Gr - Grashof number 

Gc - Solutal Grashof number 

k - Porous parameter 

Nu - Nusselt number 

Pr - Prandtl number 

Sc - Schmidt number 

R - Radiation parameter 

K - Chemical reaction parameter 

n - order of chemical reaction 
T - Temperature 
C. -Skin friction 

SS- Steady-State 

Sh - Sherwood number 

u, v - velocities in the x and y-direction respectively 

x, y - Cartesian coordinates along the plate and normal to it 

respectively 

B - Magnetic field of constant strength 
M - Magnetic field parameter 

Greek Letters 

P* - coefficient of expansion with concentration 
P - Coefficient of thermal expansion 

<t - Stefan-Boltzmann constant 

a -Scalar electrical conductivity 
p- Density of fluid 
6- dimensionless temperature 
v- Kinematic viscosity 

y - Suction parameter 

A - Variable thermal conductivity 

Subscripts 

w - condition at wall 
oo - condition at infinity 
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